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The energy, or mass scale AfsuSYi of the supersymmetry (SUSY) phase transition is, 
' as yet, unknown. If it is very high (i.e., 2> 10'^ GeV), terrestrial accelerators will not 

' be able to measure it. We determine MgusY here by combining theory with the cosmic 

microwave background (CMB) data. Starobinsky suggested an inflationary cosmological 
scenario in which inflation is driven by quantum corrections to the vacuum Einstein's 
equation. The modified Starobinsky model (MSM) is a natural extension of this. In 
the MSM, the quantum corrections are the quantum fluctuations of the supersymetric 
(SUSY) particles, whose particle content creates inflation and whose masses terminate it. 
Since the MSM is difficult to solve until the end of the infiation period, we assume here 
that an effective inflaton potential (EIP) that reproduces the time dependence of the 
cosmological scale factor of the MSM can be used to make predictions for the MSM. We 
predict the SUSY mass scale to be JV/gusY — lO-'^^GeV, thus satisfying the requirement 
that the predicted density fluctuations of the MSM be in agreement with the observed 
CMB data. 
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1. Introduction 

Cosmological models based on the anomaly-induced effective action of gravity take 
into account the vacuum quantum effects of particles in the early universe. These 
models naturally lead to inflation, which is a direct consequence of the assump- 
tion that, at high energies, all particles can be described by massless, conformally 
invariant fields with negligible interaction between them. 

In the modified version of the original Starobinsky model ^ [modified Starobin- 
sky model (MSM)], it was assumed that some of the scalar and fermion fields are 
massive. If we also assume a supersymmetric particle content, inflation is, then, 
stable during the period which starts at a sub-Planck scale, continuing almost until 
the end of inflation, when most of the sparticles decouple. Subsequently, the 
universe enters into an unstable regime, eventually undergoing a transition to the 
FRW evolution. It is to be noted that this occurs without any need for flne-tuning, 
regardless of the values of the cosmological constant A and the curvature parameter 
k. 

In fact both the cosmological constant A and the curvature k are very small 
in the context of the high-T early universe. Present observations indicate that the 
curvature fc-term is less than 10% of the energy density term in the Priedmann 
equation and that the A term is comparable to the energy density term in the 
ACDM model. Extrapolating back to the high-T early universe, the fc-term is then 
< 10~^°% and the A term < 10~^°°% of the energy density term. The fc-term is 
thus not included in our equations and the A-term is included up to Eq.(17), just 
for completeness. We set A = in the equations after Eq.(17). 

It is not possible to solve analytically the equation for the cosmological scale 
factor a{t) in the MSM. However an approximate solution for a{t) during inflation 
was obtained and conflrmed by numerical analysis with very high precision. ^ 

We assume here that an effective inflaton potential (EIP) that reproduces the 
time dependence of a{t) of the MSM during the inflationary period can be used to 
make predictions of the MSM throughout this period, including its end. We use a 
reverse engineering method, discussed in Ref. 10, to derive the EIP from a{t). 
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In § 2, wc give a brief review of anomaly- induced inflation of the MSM and the 
approximate time dependence of a(t). The effective potential and value for Msusy 
are derived in § 3. Our conclusions are presented in § 4. 



2. The approximate solution of the cosmological scale factor in the 
modified Starobinsky model 

The interesting property of the Starobinsky model is that inflation is a result of 
vacuum quantum effects of predicted (e.g., SUSY) particles. In the simplest case, 
based on the effects of massless flelds, the leading quantum phenomenon is the 
conformal anomaly. The underlying theory includes A^o scalars, N1/2 Dirac spinors, 
and A''i vectors, corresponding to the particle content of the quantum theory. They 
are not necessarily the real matter filling the universe, whose energy density is 
assumed to be negligible during inflation. The vacuum quantum effects originate 
from the virtual particles. 

Due to conformal invariance, the fields decouple from the conformal factor of 
the metric. In this case, the dominating quantum effect is the trace anomaly, which 
comes from the renormalization of the conformal invariant part of the vacuum 
action, ^^'^^ 

^vacuum — Shd + Seh , (1) 
where the first term contains higher derivatives of the metric, 

Shd = J d^xy/^ {aiC^ + a2E + aaV^R} , (2) 



and 



Seh = J d^xV^{R+2A) 



is the Einstein-Hilbert term, where ai,2,3, G, and A are the parameters of the 
vacuum action. and E are the square of the Weyl tensor and the integrand of 
the Gauss-Bonnet term, respectively: 

E = Ri.uafjR'"'"^ - 4 RafjR"^ + R^ . 
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Quantum corrections to the Einstein equation, 

R^u -^Rg^u = SttG < T^^ > - A , 

produce a non-trivial effect due to the anomalous trace of the energy-momentum 
tensor. The anomaly-induced effective action can be found explicitly. 1^-14 rpj^g 
expression for the anomalous energy-momentum tensor is 

2 ST 

T = < t; > = ff^,-— = - {wC^ + bE + cV^R) , (3) 

V 9 ^gijtv 



where w, h and c are the functions for the parameters ai, a^, as in Eq.(2), respec- 
tively: 

1 /iVo ^ ^1/2 ^ A^i\ 



6 = - 



)2 Ueo 360 180 y ' ^ ' 

(47r)2 1 180^ 30 10 J ^ ' 



and r is the quantum correction to the classical vacuum action, taking into account 
the vacuum quantum effects of the matter fields, which are free, massless, and 
conformally coupled to the metric. In terms of the new variables, and a (= 
Ina(f)), where a{t) is the cosmic scale factor and g^i, = g^v-e'^'^, the quantum 
correction to the classical vacuum action is 

f = Sc[g^,u] + j d^xy/^{waC'^ + ba{E- ^DR) + 2baAa } 

-^-^ Jd'xV^gR\ (7) 

where Sc[gp,iy] = Sclg^^v] is some unknown functional of the metric. In general, 
there is no standard method for deriving Sc[gfiv]- If we consider an isotropic and 
homogeneous metric, = g ■ [rf) , where 77 is the conformal time, the conformal 
functional Sc[gfiu\ is constant. In this case, Sc[g^ji.v\ does not depend on 0(77) and, 
therefore, does not contribute to the equations of motion [Eq.(7)] which is an exact 
one-loop quantum correction. 

The total action with quantum corrections is 

^total ~ ^vacuum ~t~ , (8) 
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which leads to the following equation of motion for a{t): 

"a d'ci a' ( 46\ aa Mp, / a a 2A\ , , 

-+3 + 5+- ---^ - + --— = 0, 9 

a a a \ c ) aa^ c \a 3/ 

where Mpi = I/VSttG = 2.44 x 10^* GeV is the reduced Planck mass. The dot 
over the variable denotes the derivative with respect to the physical time which 
is related to the conformal time r] by the relation dt = a{r])dr]. 

The last term in Eq.(9) is the trace of the standard classical Einstein equation 
when equated to (3p — p)/6c, where is the pressure (energy density) of the 
free particles and c is given by Eq.(6). There are only vacuum fluctuations (virtual 
particles) in Eq.(9) and no free particles. The other terms (the quantum correc- 
tions) become negligible if we assume that we have a FRW evolution at late times 
with a{t) ^ t^/^ in a matter dominated era, by massive particles having negligible 
pressure, for example. The equation of motion [Eq.(9)] is, then, equal to — p^/6a^c, 

2 ... 2 

instead of being homogeneous. The quantum correction terms have a , aa, a , aa , 
which are proportional to the fourth derivative of time of the cosmic scale factor a, 
while the classical Einstein term is proportional to a and a , the second derivative 
of time of a. We thus expect that the quantum correction terms will not contribute 
appreciably when the cosmic scale factor is varying slowly with time (e.g., a oc 

^2/3 

in a matter dominated universe or a oc t^^"^ in a radiation dominated universe). 
Using the FRW metric with k = 0, the solution to Eq.(9) is 

a{t) = ao ■ exp{Ht) , (10) 

where H has the form 

1/2 



H= ' . 1± Jl+(^^-^)^^ I . (11) 

v/(-327rC?)6 V ^ ^ J 

We do not consider solutions with negative H or k = ±1. The general equation of 

motion for a{t) and solutions with A; = ±1 can be found in Ref. 6. When A = 0, the 

solution Eq.(9) is 

^ = 7™i"*' '''' 

which is the solution for a{t) in the original Starobinsky model. ^ 
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The inflationary solution (10) is stable for c > and unstable for c < 0, ^ 
regardless of the cosmological constant. According to Eq.(6), the condition c > 
requires 



This enables the construction of an attractive inflationary scenario. The universe 
could start in a stable phase, such that inflation starts regardless of the initial data. 
The simplest way to provide stability in Eq.(13) is to assume that supersymmetry 
exists in the high energy region at the beginning of inflation and that it is broken 
at a lower energy since the sparticles are heavy and decouple. ^'"^ During inflation, 
H decreases due to the massive fields and, at some point, the loops of the sparti- 
cles decouple and the matter content A^o, 1/2,1 becomes modified. As a result, the 
inequality sign in Eq.(13) changes and the universe enters into an unstable inflation 
regime with an eventual transition to the FRW evolution. 

The intermediate transition epoch between inflation and post-inflation is char- 
acterized by vacuum quantum effects of both massive and massless fields. ^'^ In 
this epoch, the conformal invariance of the actions is violated by the masses and, 
therefore, we cannot use the conformal anomaly to derive quantum corrections. 
However, the conformal description of the massive theory in the framework of the 
cosmon model can be used ^^^^^ (see also Refs. 19-21 for similar analyzes and ap- 
plications of the cosmon method). A detailed discussion of the effects of massive 
fields, including quantum corrections, can be found in Ref. 6. 

A leading-log approximation for the effective action for the massive fields. 



holds in the high energy region until a "cut off" scale, defined such that at H = M* 
is reached, when a number of the sparticles decouple. The inflation then becomes 
unstable, and enters the FRW phase. 



1 < -A^i/2 + — 
3 / 18 



(13) 





(14) 
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The equation of motion for cr(t) = \na{t) in the flat case A; = 0, has the form 

b\ ...2 > .4 2AM2 

- era- - 4 - (T + 

c) c 3 c 

2 



•^• + 7aa + 4aV4(3- ^) a<7^ - 4 ^ + (X - ga - ~glA) - 



r (a + 2a') . (1 - /a) - i /V 1 = , (15) 



where 



- , 2 ^ Nfmi 



~_ 9 _ 1 Y^ iVsm^ 2 ^ A^/m) 

^ A/(87rG) 2(47r)2 ^ M^^A (47r)2 ^ M^jA " ^ ' 

In Eq.(16), the sums are taken over all fermions with a mass m/ and multiplicity 
Nf as well as over all scalars with a mass nis and multiplicity Ng. It can be seen 
that the higher- derivative terms of (14) are identical to those for the massless fields, 
as expected. ^'^^ 

Although the solution of Eq.(15) can not be performed analytically, an approx- 
imate solution can be obtained. Prom Eqs.(lO) and (11), a{t) (= lna{t)) is propor- 
tional to t, for small t. However, inflation slows down at larger t, when a{t) is no 
longer Hnear in t due to the masses of the particles. We then approximate a{t) as a 
second order term in time, a{t) = Axt + Bxt'^, where A and B are constants and 
B is negative. The equation of motion, modified by the contribution of the massive 
fields, produces the foflowing approximate solution to Eq.(15): 

where Ha = Hs/ Mpi = l/V— 2 6 is a dimensionless Hubble parameter, which de- 
pends on the particle content of Eq.(5). We use the particle content 1/2,0 = 
(12, 48, 104) of the minimum supersymmetric standard model (MSSM) since it pro- 
vides a stable infiation in the MSM and A = 0. ^ 

In a previous paper, a relatively large value for /(= 10~^) was used as a strong 
test of the accuracy of the approximate solution Eq.(17). ^ The peak value of Eq.(17) 
with this value of / was compared with the peak value of the numerical solution 
of Eq.(15) and found to differ by only one part in 10^. In Fig. 1 we plot cr(r) vs 
r(= t/\/87rtpi), where tpi = l/\/87rMpi ~ 5.3 x 10~^*sec is the Planck time, for 
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1x10^^ 2x10^^ 3x10^^ 4x10^^ 




1x10^^ 2x10^^ 3x10^^ 4x10^^ 



Fig. 1. a) Top: The dashed line shows the numerical solution of it(t s t/VSntpi), using Eq.(15), 
for A = 0, MSSM particle content and / ~ 10"" . The gray line shows the parabolic solution 
[Eq.(17)]. b) Bottom: The difference between the approximate UAprox [Eq.(17)] and the numerical 
solution (TNum [Eq.(15)]. The difference corresponds to approximately one part in 10^^. 
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/ ~ 10"^'^ [Eq.(32)] for the numerical solution [Eq.(15)] and for the approximate 
parabolic solution [Eq.(17)]. The approximate parabolic and numerical solutions 
basically coincide. The two curves differ by one part in 10^^, as seen in Fig. lb. 
Thus, for any value of / on the order of 1G~^ or less, the approximate solution is 
extremely accurate. The numerical analysis confirms the parabolic dependence of 
Eq.(17) to a very high precision up to ~ M*, where is the scale when most of 
the sparticles decouple, the inequality in Eq.(13) changes sign, and inflation becomes 
unstable. 

In § 3, we derive the reverse engineered inflaton potential that produces the ap- 
proximate solution [Eq.(17)]. In the next section, we present the theory and method 
for constructing the potential. 

3. The effective potential of the modified Starobinsky model 

Consider the dynamics of a Robertson Walker universe model with a classical scalar 
field (p{t), the inflaton and a non-interacting fluid. Following Ref. 10, we construct 
the potential V{(t)) from the set of equations 

and 

12 1 



Using the above equations, we construct the potential of the MSM: From Eq.(17), 
the Hubble parameter is 

Substituting Eq.(20) into Eq.(19), we obtain 

t{cl>) = ±^^{c^{t)-cf>o), (21) 



where |0o| > \(t>\- Choosing the positive sign in Eq.(21), we have — oo < (f) < 0. 
Substituting Eq.(21) into Eq.(20), 



i^W = ^-^(^W-0o). (22) 
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We can find V{4>) by substituting Eq.(22) into Eq.(18) or V{t) substituting Eq.(20) 
into Eq.(18): 

2" 



V{t) 



1 



(87rG) 



2(87rG) ' " i 2(87rG) 



+ 3 



t 



(23) 



In Fig. 2, we show the effective potential until the end of inflation as a function 
of T = </V87f tpi, using / ~ 10"^^ [Eq.(32)] and the MSSM particle content. The 
potential is negligible at the end of inflation {Vend/Mp^ ~ H^f). 




2x\&^ 3x10" 



4x10^^ 



Fig. 2. The potential [Eq.(23)] as a function of t(= t/VSTrtpi) for / ~ IQ-^^ and the MSSM 
particle content. 



The end of inflation occurs near the maximum of cr{t), when a{t) is no longer 
increasing exponentially with time at the mass (energy) scale M*: 



We define the dimensionless parameter 



(24) 



(25) 
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The time as a function of /x at the end of inflation is 

1 



^end — 



Haf 



JJ_ 



(26) 



Substituting ignd into Eq.(23), we have 

where Mpj = l/(87rG). Since from Eq.(32) / is smaU, V{<t)) is smaU at the end of 
inflation, compared with its initial value. 

The slow roll parameters in terms of the Hubble parameter are 



e = 2M] 



PI 



H'{cj>) 



r] = 2 Mpi 



H{<p) 

H" {<!>)- 



(27) 



H{<f>) 

The value for /i, that characterizes the end of inflation is determined by the 
condition that the slow-roll parameter of inflation is e = 1. Using this condition and 
Eqs.(22) and (27), we obtain 



The number of e-folds of inflation before t = <end is 



N 



/ H{t)dt = f7(iend) - CT{teo) 

J tp,n 



(28) 



(29) 



We are interested in A'^ ~ 60, the approximate time ^go when the observed 
density (scalar) fluctuations and the primordial gravitational (tensor) fluctuations 
from inflation were created. Substituting fond in Eq.(17), we find (Tend- From (29), 
we get aeo- Solving (17) for fgO) we obtain 

2^8^ 



*60 — 



1 



Ha 



60/ 



(30) 



Substituting tgo from Eq.(30) and pL from Eq.(28) into Eq.(23), we obtain the density 
fluctuations. 
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that is observed to be « 10~^. From the requirement that 6p/p ~ 10~^, we find 
that 

/ w 2.05 X 10-^3 . (32) 

Using this result in Eqs.(28) and (25), we obtain the scale 

« 4.5 X lO^^GeV , (33) 

which is the value of the Hubble rate -ff* when the massive particles decouple. 
Particles decouple when their masses are of the order of the temperature or of order 

i^y^Mp(^ which gives 

MsusY ^ VM^pi ~ lO^^GeV . (34) 

This value is greater than the electroweak scale (~ lO'^GeV) and consistent with 
the predicted GUT scale (~ lO" - lO^^GeV). 

4. Conclusions 

A scalar potential was constructed using the approximate solution for the time 
dependence of the cosmological scale factor of the MSM during inflation. The po- 
tential was normalized at a time 60 e-folds before the end of inflation in order 
to obtain the observed level of density fluctuations in the CMB, Sp/p 10~^. The 
mass (energy) scale of the MSM at the end of the inflation, ~ lO^^GeV, which 
we identify with the Hubble rate when the massive particles decouple, predicts a 
SUSY scale Msusy, consistent with the GUT scale Msusy ^ lO^^GeV. 
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